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OPTICALLY FAINT COUNTERPARTS TO THE /5(9-FIRBACK 170/im POPULATION: 
THE DISCOVERY OF COLD, LUMINOUS GALAXIES AT HIGH REDSHIFT 

S.C. Chapman,^ I. Smail,^ R.J. IvisoNf G. Helou,^ D.A. Dale,"* G. Lagache? 
Accepted by the Astrophysical Journal 

ABSTRACT 

We present Keck spectroscopy and UKIRT near-IR imaging observations of two 170/im- selected 
sources from the /SO-FIRBACK survey which have faint counterparts in the optical, and r — K ~ 5. 
Both sources were expected to lie at z > 1 based on their far-infrared, submillimeter and radio fluxes, 
assuming a similar spectral energy distribution to the local ultra-luminous infrared galaxy (ULIRG) 
Arp220. However, our spectroscopy indicates that the redshifts of these galaxies are z < 1\ z = 0.91 
for FNl-64 and z = 0.45 for FNl-40. While the bolometric luminosities of both galaxies are similar to 
Arp220, it appears that the dust emission in these systems has a characteristic temperature of ~30K, 
much cooler than the ~50K seen in Arp220. Neither optical spectrum shows evidence of AGN activity. 
If these galaxies are characteristic of the optically faint FIRBACK population, then evolutionary models 
of the far-infrared background must include a substantial population of cold, luminous galaxies. These 
galaxies provide an important intermediate comparison between the local luminous IR galaxies, and the 
high redshift submillimeter-selected galaxies, for which there is very little information available. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: 
starburst 



1. INTRODUCTION 

The spectral shape of the far-infrared background 
(FIRB) detected by FIRAS 100/im-4mm and DIRBE 
at 140 and 240 /im (Puget et al. 1996; Fixsen et al. 1998) 
indicates a peak at ~ 200 /im with energy comparable to 
the optical/UV background. This peak arises from opti- 
cal/UV radiation from star formation and AGN activity in 
obscured galaxies at 2: which is absorbed by dust and 
reradiated in the far-infrared. This obscured population 
of galaxies could host approximately half of the massive 
star formation activity over the history of the Universe 
(e.g. Blain et al. 1999). 

Far-infrared surveys at wavelengths close to the peak of 
the FIRB provide a powerful route for understanding the 
properties of the obscured activity in the distant Universe 
and its relevance to the formation and evolution of both 
galaxies and super-massive black holes. The FIRBACK 
(Far-InfraRed BACKground, Dole 2000) survey obtained 
wide- field imaging at 170/im with the PHOT instrument 
on-board 750 satellite in three separate regions of the sky 
chosen for low Galactic cirrus foreground. FIRBACK is 
the most reliable and deepest (cr(170/im)~ 40mJy) in- 
frared census at wavelengths near the peak in the FIRB. 
The FIRBACK sources down to 120 mJy account for about 
10% of the FIRB seen by COBE at 140-240 /im (Puget 
et al. 1999). Evolutionary models (Dole et al. 2001, La- 
gache et al. 2002) suggest that the sources identified by 
FIRBACK comprise both star-forming galaxies at low red- 
shifts, z ~ 0.1, and a population of much more luminous 
galaxies at higher redshifts. The models predict that a 
quarter of the FIRBACK sources should have z > 0.5, 



with a tail reaching beyond z = 1.5, however this is quite 
sensitive to the depth and form of the evolution assumed. 
The sources in this high-redshift tail provides the strongest 
constraints on the evolution of the population contributing 
to the peak of the FIRB at ~ 200/im. For this reason the 
identification and study of these galaxies is of particular 
interest (Sajina et al. 2002, Dennefeld et al. 2002). 

The coarse beam of ISO at 170 jam beam produces large 
uncertainties in the source positions, (100^' diameter, 99% 
error circle). However, using the empirical radio-infrared 
correlation for star forming galaxies (Helou et al. 1985), 
we can exploit deep 1.4-GHz radio data from the VLA in 
C-array configuration (Ciliegi et al. 1999) to identify the 
FIRBACK sources with positional uncertainties of only 

V' (15'' VLA beamsize). Most FIRBACK sources lying 
within the sensitive region of the VLA image are detected 
in the radio (although some are not and are therefore diffi- 
cult to identify at other wavelengths). Approximately 80% 
of the FIRBACK galaxies with radio IDs are detected in 
shallow sky surveys in the optical (DP0SS2) or near-IR 
(2MASS) (Dennefeld et al. in preparation). These 170- 
/im sources correspond to relatively low redshift galaxies. 
The remaining 20% of the radio/FIRBACK-Nl field pop- 
ulation are identified in deep UKIRT iiT-band imaging at 
17 < < 21 (Sajina et al. 2002). The faint near-IR mag- 
nitudes of these galaxies support the existence of a high- 
redshift tail in the sample and provides an efficient route 
to identifying high-redshift candidates from the 170-/im 
population. 

Additional constraints on the redshifts of the distant 
FIRBACK population are provided by submillimeter ob- 
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Fig. 1. — Multi-band images (from left to right are K, z' ^ i, r, g, U) of FNl-64 {z = 0.91 - upper panels) and FNl-40 (z = 0.45 - lower panels) 
used to identify the FIRBACK population (NET archive), and our deep UKIRT K-band imaging which reliably identifies the counterparts to 
these two sources. The data for the fainter FNl-40 has been smoothed with the PSF for visibility. The seeing at UKIRT was 0.4^^ providing 
a high resolution view of these luminous dusty galaxies. The radio source position is shown by the cross-hair and in each case is coincident 
with a very red (r — X ~ 5) galaxy component. We identify several components in each source (major components labelled Jl, J2, etc) and 
conclude that both galaxies likely show a multi-component merging structure reminiscent of local ULIRGs. Each panel is 20^^ on a side with 
North up and East to the left. 



servations of a subsample of the catalog (Scott et al. 2000; 
Sajina et al. in preparation). SCUBA photometry observa- 
tions at 850 /im, in conjunction with the 170 /im and 1.4- 
GHz fluxes, provide photometric redshift estimates (e.g. 
Carilli & Yun 2000). In this paper we present Keck 
spectroscopy and UKIRT near-IR imaging observations 
of two FIRBACK sources, FIRBACK FNl-64 and FIR- 
BACK FNl-40, which were detected by SCUBA in the 
submm (Scott et al. 2000) and whose far-IR/submm/radio 
properties indicate they potentially lie at z > 1, in the 
high-redshift tail of the FIRBACK distribution. Our 
calculations will assume a flat, A = 0.7 Universe, and 
Ho=65kms~^ Mpc~^, providing a scale of 8.4 kpc per arc- 
sec at z = 0.91 (FNl-64) and 6.2 kpc per arcsec at z = 0.45 
(FNl-40). 

2. OBSERVATIONS 

The initial identification of the FIRBACK sources 
involves searching for optical counterparts in the sec- 
ond generation Digital Palomar Observatory Sky Survey 
(DP0SS2) F-band images (close to gunn-r) of this region, 
using the accurate positions of the sources from the radio 
maps. This search highlights both FNl-40 and FNl-64 as 
optically faint (or undetected) sources, at the F < 22 limit 
of DP0SS2. 

The sources are expected to be at 2; > 1 based on photo- 
metric redshifts from the FIR/submm/radio, and within 
a degeneracy between Td and z (Blain 1999). Using the 
standard Carilli & Yun (2000) estimator for radio/850 /im 
gives 2;fni-40 = 1-1, ^fni-64 = 1-3, with a la uncer- 
tainty of 0.17. Aperture matching the observations at dif- 
ferent wavelengths (radio-15'', 850 /im-19'', 450/im-ll'', 
170 jam-lOO") could introduce errors in such redshift es- 
timates, as well as SED fitting, with confusion problems 
particularly affecting the 170 fim band. However, the very 
low source densities in the radio, 850 /im, 450 /im at the 
FIRBACK source flux levels preclude any significant con- 
tamination to the beam by interlopers with a very high 
probability. By contrast the beams are large enough to 



include all the flux from high redshift objects. To pro- 
vide complete, high resolution i^-band identifications for 
the FIRBACK sources a campaign is underway involving 
deep near-IR imaging (Chapman et al. in preparation). 

2.1. UKIRT imaging 

We observed FNl-40 and FNl-64 in the i^-band at 
UKIRT using the Fast Track Imager (UFTI) as part of the 
on-going identification campaign of FIRBACK sources. 
The source positions were taken from the radio/SCUBA 
identifications as uncovered in Scott et al. (2000), with co- 
ordinates listed in the catalog of Dole et al. (2001). Each 
source was imaged for a total of 1800 s, with individual ex- 
posures of 60 s each, reaching a limiting magnitude in a 2'' 
diameter aperture of K = 20 A (5cr) in 0.4'' seeing. Data 
were reduced using the UKIRT software pipeline oracdr 
(Bridger et al. 2000). 

The i^-band imaging of both sources is shown in Fig. 1. 
The precise mapping of the radio source catalog onto 
the optical grid suggests the astrometric error is domi- 
nated by the centroiding of the radio sources themselves 
^FWHm^ _ ^ffy rj.^^ UKIRT observations clearly 

identify near-IR counterparts at the position of the radio 
source in both cases. Indeed, we find several near-IR com- 
ponents close to the radio counterpart, on scales of < 5'' 
(<50kpc at z < 1). These components are all resolved in 
0.4'' seeing, confirming they are all galaxies. 

Based on the available optical limits we find that in both 
sources the component which is coincident with the radio 
emission, and by implication is the source of the luminous 
far-IR emission, is very red, {R — K) ^ 5 (Table 1). 

2.2. Archival INT imaging 

Archival, multi-band data is made publically avail- 
able through the Isaac Newton Groups' Wide Field 
Camera Survey Programme. We extracted U^g^r^i^z' 
band imaging for our FIRBACK galaxies, stacking mul- 
tiple exposures in the same filters with IRAF IMCOMB. 
The imaging has 0.33"/pixel, with an average 1.2" 



Chapman et al. 



3 



galaxies assuming a spectral energy distribution (SED) for 
these luminous, dusty galaxies similar to that of Arp 220. 
We now investigate this discrepency in more detail. 

3.1. Far-IR properties 

Figure 3 depicts the measured near-IR through radio 
photometry for the two FIRBACK galaxies. A template 
fit from SED library of Dale et al. (2001,2002) is overlayed, 
encompassing the mid-IR through radio wavelengths. This 
template library represents a parametrized sequence of 
SEDs (ordered in terms of their 60/100 /im restframe flux 
ratios) where a range of dust properties (temperatures and 
emissivities) are modelled to provide a composite SED, 
with increasing far-IR colors for increasing far-IR lumi- 
nosity. These SEDs provide a good representation of all 
local far-IR luminous galaxies (Dale et al. 2001,2002). The 
radio luminosity of the model SED is tied directly to the 
far-IR luminosity, using the empirical relation from Helou 
et al. (1985). 

To fit these models to our observations we simply mini- 
mize between the SED and the observations for the ra- 
dio through far-IR points, ignoring the poor signal to noise 
450 /im point, and weighting each of the points equally. In 
the case of FNl-64, the Dale et al. templates cannot simul- 
taneously fit the radio, 850/im and 170/im bands to within 
the la uncertainties, however the fit lies within the 0.2 dex 
scatter of the FIR/radio relation (Helou et al. 1985). For 
reference, we also show the best fit excluding the 170/im 
point for illustration. The excellent resulting fit with a 
slightly cooler template (lower 60/100 parameter) to the 
450 //m, 850 //m and radio points suggests that the large 
/SO/PHOT beam (100'' diameter, 99% error circle) may 
include an additional far-IR source in addition to the one 
isolated by the radio and submm. A second radio source 
with Si.4GHz < 5(7 (Ciliegi et al. 1999) lies towards the 
edge of the FIRBACK error circle, and may be contribut- 
ing to the FIR flux. 

From our fits we estimate that galaxies are relatively 
cool with rest frame 60/100 /im ratios of 0.43=b0.03 and 
0.69=b0.04 for FNl-40 and FNl-64 respectively, where flux 
errors lead to an uncertainty in the fit 60/100 /im ratio. As 
the Dale et al. SED models assume a superposition of black 
bodies of different temperatures, they can be described 
best in terms of percentage Ljr contributions from dust 
at different temperatures: Td < 15 K, < 20 K, < 31 K. 
FNl-40 has respectively 42%, 89%, 100%, while FNl-64 
has 31%, 78%, 98%. However, for direct comparison with 
other galaxies from the literature we also fit a single dust 
temperature and emissivity {/3) to the submm/far-IR data 
points and obtain: Td=25.7±0.4K and /?=1.76±0.01 for 
FNl-40; and Td=30.8zb0.7K with /3=1.68zb0.02 for FNl- 
64. These values of (3 lie within the la range of local values 
for IRAS galaxies, as fit by Dunne et al. (2000) (median 
P=1.33). We consider the effects of dust temperature fur- 
ther in subsequent sections. 

In terms of bolometric far-IR luminosity (40 to 200/im), 
using the best fitting SEDs we estimate that FNl-40 has 
a luminosity of around LpiR ~ 7 x 10^^ L©, placing it just 
in the ULIRG category. FNl-64 is more luminous, with 
LpiR 4 X 10^^ L©, placing it well into the ULIRG class 

^For completeness we note that an attempt to detect redshifted Ho; in FNl-64/ J2, using a 1800 s J-band observation with the CGS4 
spectrograph on UKIRT failed to detect a line at the relevant wavelength. 



FWHM PSF. The calibration zero points were taken from 
the archive listing for the appropriate observation date 
{www. ast. cam. ac.uk/^wfcsur/photom. html). Fig. 1 shows 
these images scaled to the UKIRT frame. FNl-64 displays 
a striking merger morphology, with components spanning 
a range in SEDs. The strong K-band Jl component ap- 
pears only as a faint tail in the optical. FNl-40 is a much 
fainter optical source, and we have smoothed the image 
with the PSF for display purposes. Some faint emission 
between Jl ane J2 in the I^z' images suggests that this 
source may also be an interaction between multiple compo- 
nents. We find that in both sources the component which 
is coincident with the radio emission, and by implication 
is the source of the luminous far-IR emission, is very red, 
{r-K)r^b (Table 1). 

2.3. Keck spectroscopy 

Having identified reliable counterparts to the 170-/im 
sources in these two fields we are then in a position to test 
the prediction that these galaxies should lie at 2: > 1. 

Spectroscopic observations of FIRBACK FNl-40 and 
FNl-64 were taken using the Echellette Spectrograph and 
Imager (ESI) on the Keck II telescope in 2001 July. In the 
echellette mode used here, ESI provides complete cover- 
age of the optical waveband from the atmospheric cut-off 
to the limit of the sensitivity of silicon CCDs: 0.32-1 /im. 
The high spectral resolution of the system, ~lA, allows for 
very good sky subtraction into the atmospheric OH forest 
at the red end of the spectrum. 

The galaxies were acquired and centered on the slit man- 
ually using the optical guide camera. In the case of FNl- 
40, we were able to center on the bright galaxy, J3 (Fig. 1), 
with the slit aligned to cover Jl, the component identified 
with the radio emission. For FNl-64, only the optically 
brighter component of the merging pair, J2, was aligned 
on the slit. 

The spectral integrations were 1800 s for each source. 
High signal-to-noise fiats and wavelength calibrations were 
taken shortly before the observations of each target, and 
fluxes were calibrated using spectra of red standard stars. 
The data were reduced with the makee software using the 
standard reduction recipe outlined in the manual (Blight- 
man 1998). The calibrated spectra for FNl-40/ Jl and 
FNl-64/ J2 (as well as the acquisition object J3) are shown 
in Fig. 2. 

The spectra of both FNl-40/ Jl and FNl-64/ J2 show 
a series strong emission lines in the red, we identify the 
strongest of these as [On] 3727, [Olll]5007 and for FNl-40: 
Yla. Based on these identifications we determine a redshift 
of ^ = 0.45 for FN1-40/J1 and z = 0.91 for FN1-64/J2.6 
In addition, we measure a redshift of z = 0.40 for the ac- 
quisition galaxy, J3, in the field of FNl-40. More details 
of the spectral line strengths from these observations are 
given in Table 1. 

3. RESULTS 

The ESI spectroscopic observations provide reliable red- 
shifts for the proposed radio-detected counterparts to the 
170-/im sources in the fields of FNl-40 and FNl-64. These 
redshifts, z = 0.45 and z = 0.91, are below those predicted 
on the basis of the far-IR/submm/radio properties of these 
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Fig. 2. — Keck/ESI spectra of FNl-40 (above) with the radio identified source Jl at z = 0.45. Inset for comparison is the optically bright 
component J3 (z = 0.40) lying on the same slit, with fiux scaled down by 100 X. While relatively close in redshift, the velocity difference 
(~13,000 km/s) is too large for them to be part of the same merging system. The Balmer ratios do not show any significant reddening. 
FN1-64/J2, the optically brighter component of an obviously merging pair, is presented below (z = 0.91). Detected lines are indicated. 
Spectral resolution is ~1A, resolving the [Oil] doublet in all cases. 



and showing that it is four times as luminous as Arp220. 
In both cases, adopting total infrared luminosities, Ltir, 
covering all the emission from 3-1100/im in the restframe, 
we estimate luminosities which are ~ 2x higher than the 
corresponding Lfir. 

We compare in Figure 4 the rest frame 60/100 ratios 
(corresponding roughly to dust temperature) and FIR lu- 
minosities for these FIRBACK galaxies with a sub-sample 
of the IRAS Bright Galaxy Sample (BGS) with 850 /im 
measurements (Dunne et al. 2000). Our extrapolation to 
60/100 is somewhat model dependent; while our observed 
170/im point measures the rest frame 100/im, our model 
fitting gives a variation in the 60/im point, dominated 
by photometric uncertainties (shown with error bars in 
Fig. 4). Fitting single dust temperature grey-bodies to lo- 
cal IRAS-selected galaxies using the 850 /im point (Dunne 
et al. 2000) has demonstrated the 60/100 ratio alone to 
slightly overestimate Td, with Td showing a tighter corre- 
lation than 60/100 with LpiR. This comparison shows that 
these two FIRBACK galaxies are colder than the major- 
ity of the IRAS BGS given their LpiR. However, the BGS 
(Seo/im >5.24 Jy) is relatively small with very few objects 
as luminous as FNl-40, FNl-64. We would like to make a 
comparison with the fainter Sgo/im 

>1.2Jy /i?^ 5 galaxies 
(Fisher et al. 1995). The local 60/100 /im distribution of 
the 1.2 Jy sample is characterized in detail in Chapman 
et al. (2002a). Placing the two FIRBACK galaxies within 
their respective far-IR luminosity class, we find that FNl- 
40 lies midway into the first quart ile of 60/100 values for 
Lfir ^ 10^^"^, while FNl-64 lies at the first quartile point 
for Lfir ~ 10^^'^. This is roughly consistent with the ex- 
trapolated relation for the BGS shown in Figure 4. We 
discuss the implications of this discovery in section 4. 
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Fig. 3. — The observed spectral energy distributions of FNl-40 
(top) and FNl-64 (bottom) from optical through radio measure- 
ments. The best fit models from the catalog of Dale et al. (2001) are 
overlaid, illustrating the colder nature of FNl-40 at z = 0.45. FNl- 
64 has hotter dust but is still relatively cool compared to Arp220 
(a ULIRG which has only 1/4 the bolometric luminosity of FNl- 
64). For Nl-64, a simultaneous fit to 170/im, 850 /im, and 1.4 GHz 
is rather poor, and an improved, cooler dust, fit to the 450 /im, 
850 /im, 1.4 GHz bands is shown as a dashed line (motivated by the 
possibility of an additional source of far-IR emission contributing to 
the large ISO-FIRBACK beam). 



3.2. Optical properties and morphologies 
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The existing optical imaging (1.2^' FWHM) reveals an 
obviously merging, r = 21.4 elongated source for FNl-64, 
and barely detected r ~ 24.5 components for FNl-40, Our 
OA" UKIRT imaging at i^-band, however, provides signif- 
icant morphological information (Fig. 1). 

The K image of FNl-64 resolves into a pair of galaxies, 
with approximately equal i^T-band flux, and the complete 
system has a K = 18.15. However, only the northern K- 
band source has a bright optical counterpart, with only a 
faint r = 23.5 tail corresponding to the southern source. 
Hence this pair comprise a blue component to the north 
and a red {r — K) ^ 5 galaxy in the south. The latter is 
coincident with the radio emission and by implication is 
probably generating the bulk of the far-IR output. The 
galaxies are separated by 2.4^' which translates into 20kpc 
at z = 0.91 and there is a bridge of emission connecting 
the two components, suggestive of tidal debris. 

FNl-40 is a faint galaxy (r = 24.5), with a possible close 
companion at 2.3'' separation (J2), and a third optically 
bright galaxy (J3) lying 3.6'' away from Jl. The radio 
centroid is coincident with component Jl, the i^-brightest 
source with {r — K) = 4.8, as labelled in Fig. 1. Our spec- 
trum identifies this as a galaxy at z = 0.45. As the J3 
source {z = 0.40) has a velocity offset of ~ 13000 km s~^ 
from Jl, it cannot be considered a part of the system. 
Both Jl and J2, however, show low surface brightness ex- 
tension (especially in the z' and i images) along the axis 
separating the two components, and we suggest that this 
probably represents a merging system. 

A comparison of the morphologies of these two FIR- 
BACK sources, multiple components separated on scales 
of - lOkpc, with local ULIRGs (Goldader et al. 2002) 
suggests a similarity with "early stage mergers" such as 
VV 114. Equally, the tendency for the radio emission (and 
therefore presumably the far-IR) to come from the red- 
der component is also shared with local ULIRG samples. 
However, there are differences: The far-IR curves for FNl- 
40/FN1-64 are fit by cool SEDs, which would normally be 
accompanied by a modest LpiR/Lopticai ratio for such early 
stage mergers (Dale et al. 2000). In contrast, the measured 
LpiR/Lopticai ratio, considering both the Jl and J2 compo- 
nents in each case, is large (20.7 and 36.2) for FNl-64 and 
FNl-40 respectively). Lopticai is calculated from the inte- 
grated fiux under a power law fit to restframe wavelengths 
0.5 /im to 1.5 /im. These values are similar to Arp220, 
a more evolved ULIRG with a hotter dust temperature 
(- 50 K). 
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Fig. 4. — The dependence of rest frame 60/100 micron ratio 
(roughly dust temperature) with LpiR for the Dunne et al. (2000) 
SCUBA observed sample of BGS IRAS sources, compared to our 
two FIRBACK galaxies, with error bars denoting the model fitting 
error from photometric uncertainties. The correlation and Icr, 2cr 
deviations are overlaid. Arp220, the ULIRG typically used for com- 
parison with high— z luminous galaxies, is identified. 

3.3. Spectral Properties and Metallicities 

Combinations of strong emission lines can be used to de- 
termine or constrain the nature of the ionizing radiation 
(Veilleux & Osterbrock 1987), the amount of reddening to 
the emission region, and the metallicity of the interstel- 
lar medium (ISM; e.g., Pagel et al. 1979; Kobulnicky et 
al. 1999). It is important to stress that these measure- 
ments only reflect the environment from which photons 
can escape and it is possible that the most active regions 
within these systems are so highly obscured that they do 
not contribute to the optical spectra used in our analysis. 
Moreover, in the case of FNl-64/ J2, our optical spectrum 
is offset from the very red object, radio emitting compo- 
nent of the system (although this has not traditionally 
limited the analysis of other similar systems, e.g. Ivison et 
al. 2001). 

The narrow lines and the resolution of the [Oil] doublet 
in both systems suggest no AGN component is present 
in either FNl-64/ J2 or FNl-40/ Jl, although without di- 
rect spectroscopic information from the far-IR dominant 
region of FNl-64/ J2, we cannot rule out an AGN origin 
within this component (Jl). The optical spectra of both 
FNl-64/ J2 and FNl-40/ Jl resemble that of luminous Hll 
regions. 

For FNl-40/ Jl, we can calculate the optical reddening 
based on the Balmer ratio, Ha/H/5. The region emitting 
Balmer lines is apparently highly reddened, equivalent to 
an Ay =3.6 mag, and consistent with the copious far-IR 
emission being centralized in the same location as the op- 
tically brighter region for which the spectra diagnostics 
have been measured. Applying this extinction to the 11/3 
inferred star formation rate (SFR) suggests a correction 
from O.7M0yr-^ to 2OM0yr-\ stih well beneath that 
extrapolated from the far-IR (16OM0yr~^). We are able 
to probe the extinction in FNl-64/ J2 through the H/?/H7, 
revealing a more modest Ay = 2.1. As stressed before this 
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Table 1 

Properties of FIRBACK FNl-64 and FNl-40 



Property 


FNl-40 


FNl-64 


A ( T9nnn^ 

XX. -TV. ^vjziuuuy 


J-U.Uc/.ZiO.U-L 


-LU.UO.^O.OO 


Dec. (J2000) 


54:28:32.6 


54:38:09.5 


redshiftspec 


0.449±0.003 


0.907±0.001 


redshiftphot 


1.1±0.2 


1.3±0.2 


U 


> 24.4/> 24.4 


23.30/> 24.6 (21.28) 


9 


> 24.7/> 24.7 


24.33/23.10 (21.91) 


r 


24.45/24.57 


23.48/22.52 (21.37) 


i 


23.96/23.81 


21.64/20.75 (19.93) 


z' 


22.31/22.40 


20.09/19.82 (18.83) 


K 


19.73/20.51 


18.63/18.71 (18.15) 


SiTOAim (mJy) 


204.7 ± 44.6 


166.2 ± 42.0 


S450Mm (mJy) 


42.0 ± 29.5 


49.7 ± 19.4 


Ssso/xm (mJy) 


6.3 ± 1.4 


5.9 ± 1.4 


b20cm (MJyj 


ooO ± oO 


OQn -L /in 


T 2 

^TIR 


6.7x10^^ 


4.0x10^^ 


1.6x10^^ 


7.7x10^^ 


[Oii]3729/3726 


1.45 


0.72 


Ha/H/3 


9.54 


n/a 


[Oiii]5007 / [iVii]6584 


1.66 


n/a 


H7/H/3 


< 0.2 


0.24 


Rl2> 


6.27 


3.56 


[Oiii]4959+5007 / [On] 


0.61 


0.95 


SFR (H/3) (reddening corrected) 


20 


416 


SFR (far-IR) 


160 


950 


FIR/Opt^ 


36.2 


20.7 


60/100 (rest frame) 


0.43 


0.69 



Quantities with divisions are quoted for components J1/J2 respectively. Magnitudes are 
calculated for 2" l2" apertures (for J1/J2). 

Quantities in brackets are for the total system in the case of the obviously merging FNl-64. 
3cr photometry limits are (t/, g, r, i, 2' , K) = (24.4,24.7,24.6,23.9,22.4,20.5). 

Sfir = (1. 26x10"^^) * (2.58xS6o + Sioo) [W/m^]; Lfir computed using a flat, A = 0.7 cosmology. 
Ltir covers 3)Ltm to llOOyLtm as defined in Dale et al. (2001a). 

The obscuration ratio between FIR luminosity and optical/near-IR luminosity, as described in the text. 



constraint is not for the major source of far-IR emission in 
this system (Jl). However, the H/? imphed SFR from J2 
(416M0yr~^), is comparable to the far-IR derived SFR 
from Jl (95OM0yr-i). The FNl-64 system is therefore 
reminiscent of other high— z ULIRGs where the optically 
bright component provides an insight into the luminosity 
of the obscured portion (e.g. Ivison et al. 2001 - see dis- 
cussion for further examples). 

To measure the ISM metallicity, a rigorous determina- 
tion requires knowledge of electron temperature derivable 
only from intrinsically weak lines. However, metallicity di- 
agnostic line ratios based on stronger lines have been em- 
pirically calibrated. In particular, the i?23 parameter, i^23 
= ([Oii]3727 + [Oiii]4959+5007)/H/5 (Pagel et al. 1979), 
has been calibrated against metallicity, with an intrinsic 
scatter of only 0.2 dex. It is a weak function of the ion- 
ization ratio [Oiii]4959+5007/[Oii]3727. A reversal in i?23 
occurs at Z ~ O.SZ© due to cooling effects, so a low- and 
a high- metallicity solution are associated with most values 
of i?23- This degeneracy can however be broken using the 
[Olll]5007/[Nll]6584 ratio (Kobulnicky et al. 1999). 

The metallicity can be used as an indicator of a galaxy's 
evolutionary state, and as a constraint on possible present- 
day descendants. Carollo & Lilly (2001) have investigated 
the ISM metallicities of H/5-luminous field galaxies in the 
0.5 < 2; < 1.0 redshift interval from the Canada- France 
Redshift Survey (CFRS). These galaxies fall in the i?23 
versus [Olll]/[Oll] plane in locations that are occupied by 
galaxies in the local Jansen et al. (2000) sample. At both 
epochs, galaxies selected to have the same H/3 luminosities 
exhibit the same range of i^23 and [Olll]/[Oll]. We find 
that the galaxies FN1-40/J1 and FN1-64/J2 are within 



our la error of the median properties of the CFRS sources 
in the [Olll]4959+5007/[Oll]3727 versus R23 plane. For 
the lower redshift FNl-40 source, we can use the detected 
[Nll6584] to break the metal degeneracy of the R23 indica- 
tor. The [Om] 5007/ [A/"//] 6584 < 2 necessitates metallic- 
ities Z > O.5Z0, placing the source on the upper branch 
of the relation, with a metallicity close to the solar value. 

While the Carollo & Lilly (2001) galaxies in the 0.5 < 
2: < 1.0 range are bright star forming galaxies, they fall 
an order of magnitude below the bolometric luminosity 
of our two FIRBACK sources under study (several Car- 
ollo & Lilly 2001 galaxies are detected in the CFRS 15/xm 
ISO sample of Flores et al. 1999). While it is therefore 
interesting to compare R23 metallicities for these objects, 
we stress that the R23 indicator has only been calibrated 
for low obscuration regions, and may not be a valid surro- 
gate for metallicity in the highly obscured environments of 
ULIRGs. In particular, we are only measuring i^23 from 
regions which happen to be visible from a system which is 
likely to be generally obscured. 

4. DISCUSSION 

At z < 1, sources observed at both 170 and 850 fim will 
exhibit similar characteristics as a function of dust tem- 
perature (Blain 1999); the coldest and most dusty galaxies 
will have the greatest fiux densities for an equivalent far- 
IR luminosity. Beyond 2: ~ 1, the 170/im band lies blue- 
ward of the peak in the restframe black body and dust 
temperatures T^^ < 100 K no longer have significant ef- 
fects on the observed source fiux. A fiux-limited survey 
such as FIRBACK is thus expected to be biased towards 
cooler galaxies at higher redshifts. While we do not yet 
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have the statistics to quantify the actual selection rate, we 
note that the identification of the unusually cold ULIRGs 
FNl-40 and FNl-64 are in agreement with this prediction: 
we find cooler galaxies than expected for the ~ 10^^ 
bolometric luminosities. 

Models of the FIRBACK population (Dole et al. 2001) 
have evolved the local luminosity function by boosting only 
the ULIRG contribution with redshift, assuming an aver- 
age SED template for ULIRGs. This scenario is able to fit 
the FIRB and the FIRBACK counts, inferring that 30% of 
these sources lie at 0.5 < z < 1.0, and 10% at 1 < z < 2.5 
(Dole et al. 2001). However, the assumption that an aver- 
age ULIRG template can represent the FIRBACK source 
population is flawed as it does not take into account the 
much colder sources which will be preferentially selected at 
170/im for z < 1. Other recent models of the far-IR galaxy 
population (Franceschini et al. 2001; Chary & Elbaz 2001) 
do not incorporate a cold luminous galaxy template and 
cannot directly acount for numerous sources such as FNl- 
40/FN1-64. Revisions to the Dole et al. model (Lagache 
et al. 2002) include a population of colder luminous ob- 
jects based on the SED shape of the less luminous local 
FIRBACK sources. 

The actual redshifts of both FN1-40/FN1-64 {z = 
0.45/z = 0.91) lie significantly below the z = 1.1/ z = 1.3 
predicted from the submm and radio, using the Carilli & 
Yun (2000) estimator (with la uncertainty of 0.17) which 
is also derived from the average ULIRG SED. Moreover, if 
the majority of FIRBACK sources were to span the range 
in properties observed in these two galaxies, the average 
properties deduced from an Arp220 template would be 
severely skewed to higher redshifts. We chose these two 
sources for Keck follow-up at random from the FIRBACK 
sources where the sub- mm/radio ratio suggested a high 
{z > 1) redshift. As a consequence, the distinct possibility 
exists that almost none of the radio identified FIRBACK 
sources may lie at z > 1, dependent on the mid- and high- 
redshift distributions in dust temperature. 

The measurement of the 850/1.4 ratio for both these 
sources does allow a constraint on its use as a redshift 
indicator. Clearly the ratio for these sources must lie sig- 
nificantly above the canonical value for ULIRGs found by 
Carilli & Yun (2000), as shown in Fig. 5. The existence of 
such cooler luminous sources may affect seriously the inter- 
pretation of the high redshift SCUBA selected population 
(see also Eales et al. 2000). 

Cold yet luminous dusty sources suggest large masses 
of dust heated at moderate intensity, rather than small 
amounts of dust in a compact configuration subjected to 
extremely intense radiation fields. This latter scenario is 
clearly supported by the spatial and SED properties of 
the nearby ULIRG Arp220 (Soifer et al. 1984). While one 
could invoke evolving or novel dust properties to explain 
cooler high-luminosity sources like FNl-40 and FNl-64, 
there is no evidence to support such unusual dust proper- 
ties, especially that the SEDs of both FNl-40 and FNl- 
64 are fitted quite well with the Dale et al. models from 
the local Universe. For objects like FNl-40 and FNl-64 
however, the data suggest a spatially distributed star for- 
mation activity heating a more extended ISM. Arp 302 
represents a nearby system with a 60/100 color as cold as 
FNl-40 and Lfir=5x 10^^ L©. Lo, Gao & Gruendl (1997) 
have shown that luminous CO emission traces the entire 



optically defined disks of both merging components, cover- 
ing 23kpc and 11 kpc respectively, and we might imagine a 
similar configuration in FNl-40. The main difficulty then 
is to explain the high extinctions suggested by the large 
infrared-to-visible ratios in these systems. Since we know 
little of their detailed morphology or geometry, we cannot 
exclude the possibility that the initial conditions and the 
details of the encounters (strong tidal interactions rather 
than direct collisions or mergers) result in obscured yet 
distributed rather than strongly concentrated star forma- 
tion activity. 




1 2 3 4 5 



Redshift 

Fig. 5. — The 850/im/1.4 GHz ratio as a function of redshift 
for galaxies of various activity levels, ranging from normal spi- 
rals to starbursts, with the Carilli &: Yun (2000) relation shown 
for ULIRGs. Our two cold FIRBACK sources are shown double 
stars. Other star symbols depict sub-mm sources with no sign 
of AGN activity (in increasing redshift order: HRIO-Dey et al. 
1999, SMM J14011+0252-Ivison et al. 2001, MMDll-Chapman et 
al. 2002b, SA22-blobl-Chapman et al. 2001). For comparison we 
plot 4 sub-mm selected sources identified as AGN with filled circles 
(Small et al. 2002) and 4 BAL quasars which were sub-mm detected 
(double circles - Lewis & Chapman 2002). 

4.1. Comparison with higher redshift SCUBA populations 

Figure || shows how these FIRBACK sources com- 
pare to other submm luminous sources with measured 
spectroscopic redshifts. Note that while the FIRBACK 
sources appear to be pure starburst, there are very few 
spectroscopically identified sub-mm galaxies at high red- 
shift for comparison which do not show AGN signatures 
(SMM J14011+0252 at z = 2.56 - Ivison et al. 2001, West- 
phal MMD 11 at z = 2.99 - Chapman et al. 2002b, HRIO 
at ^ = 1.44 - Dey et al. 1999). Both FNl-64 and FNl- 
40 appear morphologically similar to the SCUBA galaxies 
SMMJ14011+0252 and Westphal MMD 11, having a red 
{R — K) source identified to the far-IR emission along with 
a bluer companion. This is similar to several local ULIRGs 
(e.g. 08279+1372 - Surace & Sanders 2000; Goldader et 
al. 2002), and may suggest that all of these galaxies are 
generating their far-IR output by similar mechanisms. 

As these two sources lie midway between the local IRAS 
population, and the high redshift SCUBA-selected popula- 
tion (e.g. Small et al. 2002), understanding their properties 
provides a crucial benchmark for continued progress on the 
evolution of IR-luminous galaxies. FNl-40 and FNl-64 
are representative of the optically faint FIRBACK sample 
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with radio and SCUBA detections (~30% - Scott et al. 

2000) , comprising -8% of the radio identified FIRBACK 
survey. In terms of the bright (Ssso/um > 6mJy), blank 
field submm sources (e.g. Small et al. 2001, Scott et al. 

2001) in this area on the sky, they would represent only a 
small fraction of the 850 //m counts 3%). 

The strong bias to colder dust temperatures in sub-mm 
and far-IR selected samples occurs because the selection 
bands fall redward of the grey-body peak. At an observed 
170 /im, this implies the cold bias out to z ~ 1 discussed 
earlier. However, for 850 /im, this bias is inherent in the 
entire plausible range of galaxy evolution, out to z ^ S 
(e.g. Blain 1999). In this sense, the wide field SCUBA 
surveys will preferentially uncover the coldest sources for 
a given luminosity. 

We would like to use the existence of these cold, lumi- 
nous sources to place a constraint on the dust temperature 
distribution at high redshifts. If the important driver of 
the dust temperature is the radiation field then we must 
compare at similar luminosity systems. Figure 4 suggests 
that these luminous FIRBACK sources are within the 
coldest 16% of the local population if we cut at the same 
luminosity. If we propose that this local population of IR- 
luminous objects undergoes a strong luminosity evolution 
of the form (1 + z)^ out to z = 1.5, we can determine the 
relative numbers of cold and warm objects isolated by our 
selection function. This form of evolution has been shown 
to match the observe space density of high— z ULIRGs to 
the local population (e.g. Blain et al. 1999; Chapman et 
al. 2002c). 

We start with the FIRBACK survey limit, S170 > 
120 mJy. We add the requirement that S850/S1.4 give 
^iy^ > 2, the CariUi & Yun (2000) criterion we em- 
ployed to select high redshift objects for spectroscopic fol- 
lowup with Keck. This selection function implies that we 
begin to detect objects with Td < 30K at z > 0.35, 
and Td < 50K at ^ > 1. We find 35% more FIR- 
BACK detectable objects with Td < 30K than with 
30 K< Td < 50K in the accessible volume from z = 0.35 
to z = 1.6, where we can no longer detect a 50K source. 
As our two detected sources have Td=26K and 32 K, our 
35% bias towards selecting Td < 30 K sources implies that 
the outcome of our experiment deviates by less than la 
from the expected outcome. 

We can conclude that the dust temperature distribution 
at higher redshifts is consistent with that observed locally 
(or alternatively, that the characteristic temperature of a 
ULIRG doesn't increase dramatically out to z ^ 1). Note 
that this is not necessarily the expected result. For in- 
stance, different temperatures distributions for similar lu- 
minosity ULIRG populations at 2: = and z = 1 could 
arise due to different geometry of the active region, or 
different chemical properties of the dust. Complete red- 
shift distributions for the FIRBACK population will be 
required to carefully test the dust temperature distribu- 
tions at low and high redshift. A more detailed model, 
taking acount of the local bivariate luminosity function 
in Ltir and Seo/im/Sioo/im, studies directly the effect of 
the complete distribution of dust properties on the high 
redshift far-IR population (Chapman et al. 2002a). 



5. CONCLUSIONS 

We have obtained Keck spectra and UKIRT high spatial 
resolution near-IR imagery for two of the proposed high- 
est redshift sources from the FIRBACK-Nl 170/im survey. 
This survey is currently the most sensitive probe of the 
properties of dusty galaxies at the peak of the FIRB. 

We find that the redshifts of counterparts to the 170 fim 
sources confirm that both sources are ULIRGs, but that 
their redshifts are significantly lower than implied by fit- 
ting a typical ULIRG SED to their far-IR/submm/radio 
SEDs. This indicates that they have cooler dust tem- 
peratures, T^^ ~ 30 K, than the canonical ULIRG values 
{Td ~ 50 K). The sources both show morphologies sugges- 
tive of early stage mergers, similar to recently identified 
SCUBA galaxies, and also many local ULIRGs. 

For FNl-40, we are able to examine optical line diagnos- 
tics in the region of far-IR emission. We measure a large, 
Av ~ 4 mag, extinction, indicating a large correction to 
the H/3 estimated SFR. The metallicity in this region, as 
probed by the R23 and indicator and OIII/NII ratio, is 
found to be typical for starforming galaxies 2 orders of 
magnitude less luminous (close to solar). 

The detection of these galaxies in the submm by SCUBA 
allows a constraint on the 850/1.4 redshift indicator. We 
find a relation for these galaxies lying significantly above 
the relation found by Carilli & Yun (2000), due to their 
cooler dust temperatures. The existence of such cooler, 
but still luminous sources may thus affect the interpreta- 
tion of the SCUBA selected population, if they are numer- 
ous at high redshift. Recent models of the far-IR galaxy 
counts, relying on the rapid evolution of a hot ULIRG 
component to the local luminosity function, may predict 
an erroneous redshift distribution without taking into ac- 
count sources such as FNl-40 and FNl-64. 

We have turned the discovery of these two cold, lumi- 
nous galaxies into an estimate for the dust temperature 
distribution at 2: ^ 1 for the ULIRG population. We 
conclude that our 170/im, 850 /im, and 1.4 GHz selection 
function isolates approximately 35% more cold sources 
(< 30 K) than warmer sources (30 <Td < 50) over the 
detectable redshift range, under the assumption that the 
temperature distribution at 2: = and 2: = 1 is similar. 
The relative proportion of the cold population at low/high- 
z can therefore be understood as consistent with the local 
60 /im/100 /im distribution. 

Thanks to A. Blain for stimulating conversations about 
this source population. We acknowledge the useful com- 
ments of an anonymous referee, which helped improve the 
text. The INT WFC data was made publically available 
through the Isaac Newton Groups' Wide Field Camera 
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ated on the island of La Palma by the Isaac Newton Group 
in the Spanish Observatorio del Roque de los Muchachos 
of the Instituto de Astrofisica de Canarias. UKIRT is op- 
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